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The lattice QCD simulation of the quark-gluon mixed condensate 
g{qa^TjG fj,Tjq) at finite temperature and the phase transition of QCD 
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The thermal effects on the quark-gluon mixed condensate g{qa^vG^vq) , which is another chiral order parameter, 
are studied using the SU(3)c lattice QCD with the Kogut-Susskind fermion at the quenched level. We perform 
the accurate measurement of (;((7(j^^G^^g) as well as (gg) for Q ^ 500MeV. We observe the sharp decrease 
of both the condensates around Tc — 280MeV, while the thermal effects below Tc are found to be weak. We also 
find that the ratio mo = g {qa ^,yG ^^q) / {qq) is almost independent of the temperature even in the very vicinity of 
Tc, which indicates that the two condensates have nontrivial similarity in the chiral behaviors. We also present 
the correlation between the condensates and the Polyakov loop to understand the vacuum structure of QCD. 



In order to understand the nature of QCD and 
the relation to hadron phenomenology, the non- 
perturbative nature is one of the most important 
characteristics of QCD. The nonperturbative phe- 
nomena, such as spontaneous chiral-symmetry 
breaking and color confinement, also bring the 
rich structure to the QCD phase diagram. For 
example, at high temperature, QCD is believed 
to exhibit phase transition into QGP, where chi- 
ral symmetry is restored and the color is decon- 
fined. To realize these phenomena, the RHIC ex- 
periments are in progress, which attempt to pro- 
duce QGP in the laboratory. For the theoreti- 
cal study of the finite temperature QCD and the 
phase transition of QCD, we focus on the thermal 
effects on the condensates. In fact, condensates 
directly characterize the nontrivial QCD vacuum 
and thus can indicate the change of the vacuum 
structure at finite temperature. 

Among various condensates, we study the 
quark-gluon mixed condensate g{q<J^vG^^q) at fi- 
nite temperature as a relevant physical quantity 
on the chiral structure of the QCD vacuum. Here, 
the mixed condensate is another chiral order pa- 
rameter of the QCD vacuum, since the chirality of 
the quark in g{qcF^uG^yq) flips as g{qcr^^G^i,q) = 
g{qR{o-)_c,yG^^)qL)+g{qLicr^^G^^)qR). We empha- 
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size that g{qa^yGf^uq) characterizes different as- 
pect of the QCD vacuum from {qq), because 
g{qo^vGfj,uq) reflects the direct correlation be- 
tween color-octet components of q-q pairs and the 
spontaneously generated gluon field, while {qq) 
reflects only the color-singlet q-q components. 

We further note that the mixed condensate 
is important quantity in hadron phenomenology 
through the QCD sum rule framework. In fact, 
it is known that g{qafiuG^uq) has large effects in 
the QCD sum rule especially for baryons, such 
as A^-A splitting P and parity splitting Re- 
cently, it is also shown P| that g{s(T^vGfius) / (ss) 
is a key quantity for the prediction on the parity 
of the penta-quark baryon, 9+ (1540) 

To study the thermal effects on g{qa^i,G^vq}, 
we use lattice QCD Monte Carlo simulation, 
which is the direct and nonperturbative calcula- 
tion from QCD. So far, the lattice QCD stud- 
ies for g{q(Jfi^G^yq) have been performed at zero 
temperature by a pioneering but rather prelim- 
inary work Uj, and by the recent works of our 
group jBj using the KS fermion, and of another 
group |3| using the Domain- Wall fermion. At 
finite temperature, however, there has been no 
result on g{q(J^^G^vq) except for our early re- 
ports |S|. Therefore, we present the extensive re- 
sults of the thermal effects on g{qa^^G^vq) with 
the analysis near the critical temperature 9 . 

We evaluate the condensates g{qcf^i,jGi^ivq) as 
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well as (qq) using the SU(3)c lattice QCD at the 
quenched level. The Monte Carlo simulations are 
performed with the standard Wilson action for 
P = 6.0,6.1 and 6.2. The lattice units are ob- 
tained as fl-i = 1.9, 2.3 and 2.7GeV for /3 = 
6.0, 6.1 and 6.2, respectively, so as to reproduce 
the string tension ^ = 427MeV [5]. We per- 
form the calculation at various temperatures as 
0<,T <, SOOMeV using the following lattices, 

(i) (3 = 6.0, 16^ X Nt [Nt = 16, 12, 10, 8, 6, 4), 

(ii) (3 6.1, 20^ X Nt [Nt = 20,12,10,8,6), 

(iii) /3 = 6.2, 243 X Nt {Nt = 24, 16, 12, 10, 8). 
We generate 100 gauge configurations for each 
lattice. Moreover, we generate 1000 gauge con- 
figurations in the vicinity of the phase transition 
point, namely, 20^ x 8 at /3 = 6.1 and 24^ x 10 
a.t f3 — 6.2, because the fluctuations of the con- 
densates get larger near T^. For the lattices at 
T > Tc, we only use the gauge configurations 
which are continuously connected to the trivial 
vacuum — 1. For the fermion action, we em- 
ploy the KS fermion to preserve the explicit chiral 
symmetry for the quark mass m — Q, which is 
desirable for the study of chiral order parameters. 

In the calculation of the condensates on 
the lattice, we use the formula as a^{qq) = 
-3E/Tr[((7/(a;)g/(a;))], a'^giqa^.G^^.q) = 
-lEf, Tr [(g/(a;)q-/(x))a,,.G;f*(x)] , where 
SU(4)/ quark-spinor fields, q and q, are con- 
verted into spinless Grassmann KS fields and 
the gauge-link variable in the actual calculations. 
The more detailed formula are given in Ref. 6_. 

In each configuration, we measure the conden- 
sates on 16 {f3 = 6.0) or 2 (/3 = 6.1, 6.2) different 
points which are taken so as to be equally spaced 
on the 4-dimensional lattice |6I9| . Therefore, we 
achieve high statistics as 1600 data at /3 = 6.0 
and 200 data at /3 = 6.1 and 6.2. Note that for 
the lattices of 20^ x 8 (/3 = 6.1) and 24^ x 10 
(/3 — 6.2), we obtain 2000 data, which guarantees 
reliability of the results even in the vicinity of Tc. 

We calculate the condensates at each quark 
mass of m ~ 20, 35, 50 MeV. At each tempera- 
ture, we observe that both the condensates show 
a clear linear behavior against m, and therefore 
we fit the data with a linear function and deter- 
mine the condensates in the chiral limit. We es- 
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Figure 1. g{qu^uG^vq)T normalized by 

g{qo'^uGf_ii,q)T=o plotted against temperature T. 
The vertical dashed line denotes the critical tem- 
perature 
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Figure 2. ml{T) = g{q(ytj.uG^iyq)T / {qq)T nor- 
malized by mg(T — 0) plotted against T. This 
result indicates the same chiral behavior between 
g{q<Jf,i,Gf,„q)T and {qq)T- 



timate the statistical error to be 5-7% level using 
jackknife error method. The finite-volume arti- 
fact is estimated to be about 1% level, through a 
check of boundary condition effects |6I9| . 

We evaluate the thermal effects on each con- 
densate, g{q(J ^^G ^v<l) or {qq), by taking the ratio 
between the values at finite and zero tempera- 
tures. Note that the renormalization constants 
cancel in these ratios. In figure we plot the 
thermal effects on g{q(Jf^i,G^j,vq) . We find a drastic 
change of g{q<J^vG^^q) around the critical tem- 
perature Tc ~ 280MeV. We obtain Tc/y^ = 
0.64(4), which is consistent with the coincidence 
of the confinement/deconfinement phase transi- 
tion and chiral-symmetry restoration. We also 
find that the thermal effects on g{q(J^,yGfj_^q) are 
remarkably weak for T <, 0.9 Tc. The same non- 
trivial features are also found for {qq). 

We then quantitatively compare the thermal 
effects of g{q(JfivG^^q) and {qq) by plotting the 
ratio ■ml{T) = g{qaf_,„Gf_t^q)T/{qq)T against T. 
In figure [3 we observe that TO§(r) is almost in- 
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Figures. Time history of g{qa^,^G^^q)T (upper) 
and the Polyakov loop (lower) for the 500 config- 
urations on 24^^ X 10 {(3 — 6.2) lattice near Tc- 
This result indicates strong correlation between 
these two quantities. 

dependent of the temperature, even in the very 
vicinity of Tc- This nontrivial result can be inter- 
preted that g{qa^vG^i,q)T and {qq)T have simi- 
larity in the chiral behaviors. 

To understand this nontrivial similarity, we ex- 
press the condensates as {qq) = y J d\' , 

g{q(^i,uG^,uq) = f / d\' ^2+^2 {>^'W^uG^,t,\\'), 
where |A) denotes the eigenvector of the Dirac op- 
erator and p{\) the spectral density on A JO]. Us- 
ing these formula, the common thermal behavior 
is understood in the way that the thermal effects 
are dominated by p(A)|a=o and (A|(Tp^Gp^|A)|A=o 
has remarkably weak dependence on T. 

This result further implies that the thermal 
effects hardly appear in the local gluon field 
strength, but appear in the global structure of 
the vacuum such as topological quantities. To 
proceed, we plot the time history of g{qafj^uG^j,vq) 
and the Polyakov loop in figure 13 and find the 
strong correlation between g{q<Tf^^G^j,^q) and the 
Polyakov loop. For the gauge configuration with a 
small (large) value oi g{q(j^j,^G^j,^q), the Polyakov 
loop takes a large (small) value. One of the pos- 
sible explanation is that the thermal effects are 
dominated by the quark propagation making a 



circuit in time direction, which can directly re- 
flect the topological structure of the vacuum. 

In summary, we have studied thermal effects on 
gigf^fiiyG^i^q) using the SU(3)c lattice QCD with 
the KS fermion at the quenched level. We have 
observed a clear signal of chiral restoration as a 
sharp decrease of g{qu^vG^^q) as well as (gq), 
while the thermal effects have been found to be 
small for T <, 0.9 Tc. We have also found that 
mg(T) = g{qa^i,G^uq)T/ {qq)T is almost indepen- 
dent of T in the entire region up to Tc, which indi- 
cates that these chiral condensates show a com- 
mon thermal behavior. The strong correlation 
between the condensates and the Polyakov loop 
has been also observed. For further studies, full 
QCD lattice calculations are interesting to ana- 
lyze dynamical quark effects on the condensates. 

REFERENCES 

1. H.G. Dosch, M. Jamin, and S. Narison, Phys. 
Lett. B220 (1989) 251; W-Y.P. Hwang and 
K.-C. Yang, Phys. Rev. D49 (1994) 460. 

2. D. Jido, N. Kodama and M. Oka, Phys. 
Rev. D54 (1996) 4532; D. Jido and M. Oka, 
"hep-ph/9611322 

3. J. Sugiyama, T. Doi and M. Oka, Phys. Lett. 
B581 (2004) 167. 

4. T. Nakano et al., Phys. Rev. Lett. 91 (2003) 
012002. 

5. M. Kremer and G. Schierholz, Phys. Lett. 
B194 (1987) 283 . 

6. T. Doi, N. Ishu, M. Oka and H. Suganuma, 
Phys. Rev. D67 (2003) 054504; Proc. of 
"Quark Confinement and the Hadron Spec- 
trum V" , (World Scientific, 2003) p.381 

7. T.W. Chiu and T.H. Hsieh, Nucl. Phys. B673 
(2003) 217. 

8. T. Doi, N. Ishii, M. Oka, H. Suganuma, Nucl. 
Phys. A721 (2003) 934; Prog. Theor. Phys. 
Suppl. 151 (2003) 161; Nucl. Phys. B129 
(Proc. Suppl.) (2004) 566; Proc. of "Color 
Confinement and Hadrons in Quantum Chro- 
modynamics" , (World Scientific, 2004) p.398. 

9. T. Doi, N. Ishii, M. Oka and H. Suganuma, 
Phys. Rev. D70 (2004) 034510. 

10. T. Banks and A. Casher, Nucl. Phys. B169 
(1980) 103; S. Hands, J.B. Kogut and A. Ko- 
cic, Nucl. Phys. B357 (1991) 467. 



